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>• ' Abstract 

en . Recent T2K results indicate a sizeable reactor angle ^13 which would rule out exact 

tri-bimaximal lepton mixing. We study the vacuum alignment of the Altarelli- 
Feruglio A4 family symmetry model including additional flavons in the 1' and 1" 

^—N . representations and show that it leads to trimaximal mixing in which the second 

column of the lepton mixing matrix consists of the column vector (1, 1, 1) /v3, with 
a potentially large reactor angle. In order to limit the reactor angle and control the 
higher order corrections, we propose a renormalisable 5*4 model in which the 1' and 

K> , 1" flavons of A4 are unified into a doublet of S4 which is spontaneously broken to A4 

^1 by a flavon which enters the neutrino sector at higher order. We study the vacuum 

alignment in the 5*4 model and show that it predicts accurate trimaximal mixing 
with approximate tri-bimaximal mixing, leading to a new mixing sum rule testable 
in future neutrino experiments. Both A4 and 5*4 models preserve form dominance 
and hence predict zero leptogenesis, up to renormalisation group corrections. 
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1 Introduction 

Recently T2K have published evidence for a large non-zero reactor angle |T] which, when 
combined with data from MINOS and other experiments in a global fit yields [2], 

ei3 = 8°±1.5°, (1.1) 

where the errors indicate the one a range, although the statistical significance of a non- 
zero reactor angle is about 3cr. 

If confirmed this would rule out the hypothesis of exact tri-bimaximal (TB) mixing [3]. 
However other schemes such as trimaximal (TM) mixing remain viable 



[/tm = P' -^cos'i^-^sin^e-*^ ^ ^ cos ^ - ^ sin ^ e*'' P, (1.2) 

\ — 4= cos?? + 4= sint^e"*^ 4= —4= cos?? — ^sin-i^e*''; 

where -j=sm'd = sin 6*13, P' is a diagonal phase matrix required to put [/tm into the 

PDG convention [9], and P = diag(l, e*"2", e* 2 ) contains the usual Majorana phases. In 
particular TM mixing approximately preserves the successful TB mixing for the solar 
mixing angle 6*12 ~ 35° as the correction due to a non-zero but relatively small reactor 
angle is of second order. Although TM mixing reduces to TB mixing in the limit that 
f? ^ 0, it is worth emphasising that in general TM mixing involves a reactor angle ^13 
which could in principle be large or even maximal (e.g. 45°). The observed smallness of 
the reactor angle 6is compared to the atmospheric angle 623 ~ 45° and the solar angle 
612 ~ 34° [2] is therefore not explained by the TM hypothesis alone. Clearly the relative 
smallness of the reactor angle can only be explained with additional model dependent 
input. We shall show that, although an A4 family symmetry can provide an explanation 
of trimaximal mixing, the observed smallness of the reactor angle as compared to the 
atmospheric and solar angles may be naturally explained by an 5*4 family symmetry. 

In the original A4 models of TB mixing Higgs fields [IQj or flavons pTl[T2] transforming 
under A^ as 3 and 1 but not 1' or 1" were used to break the family symmetry and to lead 
to TB mixing. However there is no good reason not to include flavons transforming as 1' 
or 1" and once included they will lead to deviations from TB mixing [T31[I1] in particular 
it has been noted that they lead to TM mixing [15]. However, as remarked above, TM 
mixing by itself does not account for the smallness of the reactor angle, and in particular 
the model in [15] does not provide any explanation for this. Another aspect is that higher 
order operators in A4 models lead to deviations from TB mixing in a rather uncontrolled 
way [12], and in such models any simple TM structure may be washed out unless the 
theory is promoted to a renormalisable one in which the higher order operators are under 
control. For example, this was done for two A^ models of TB mixing in [16], and one may 
follow an analogous strategy including also 1' or 1" flavons. 

In this paper we make a first study of the vacuum alignment of an ^4 family symmetry 
model of leptons with additional flavons in the 1' and/or 1" representations and show that 
it leads to accurate trimaximal mixing at leading order (LO). In order to constrain the 
reactor angle and control the higher order corrections, we then propose a renormalisable 
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5*4 model of leptons in which the 1' and 1" flavons of ^4 are unified into a doublet of 
5*4 which is spontaneously broken to A4 by a flavon which enters the neutrino sector at 
higher order. We study the vacuum alignment of the 5*4 model and show that it predicts 
accurate trimaximal mixing with approximate tri-bimaximal mixing, leading to a new 
mixing sum rule testable in future neutrino experiments. We also remark that both y44 
and 5*4 models preserve form dominance [17] and hence predict zero leptogenesis [18], up 
to renormalisation group corrections. 

The layout of the rest of the paper is as follows. In Section [2] we first revisit the 
Altarelli-Feruglio A4 model [12] with regard to the possibility of generating deviations 
from TB mixing which respect TM mixing. As the obvious ideas do not yield TM mixing, 
we then consider the model augmented by extra fiavons in the 1' and/or 1" representation. 
We study the vacuum alignment and show that this model leads to TM mixing, with an 
unconstrained reactor angle. In Section [3] we propose a renormalisable 6*4 model of leptons 
and study its vacuum alignment leading to accurate trimaximal mixing with approximate 
tri-bimaximal mixing. In Section H] we give an analytic discussion of the perturbative 
deviations to TB mixing arising from any TM model with a physical reactor angle leading 
to a new mixing sum rule testable in future neutrino experiments. Section E] concludes the 
paper. Appendix A contains Clebsch-Gordon coefficients for S4 and A4, and Appendix B 
describes a perturbative diagonalisation of the right-handed neutrino mass matrix. 

2 A4 models of trimaximal mixing 

2.1 The Altarelli-Feruglio A4 model of tri-bimaximal mixing 

The original and well studied Altarelli-Feruglio (AF) model of lepton masses and mixings 
[12] is formulated as an effective theory, defined purely by the particle content and the 
symmetries. There exist two versions, one with right-handed neutrinos and one without. 
For definiteness we will only consider the former which makes use of the elegant seesaw 
mechanism to obtain effective light neutrino masses around the eV scale. The particle 
content and the symmetries of the AF model we consider is presented in Table [1], including 
the ^ flavon, but excluding the ^' and ^" flavons which we shall consider later. In addition 
the vacuum alignment in the AF model requires a further auxiliary flavon singlet ^ which 
does not acquire a VEV and is not shown in the table. Particles with a U{1)r charge of 
2 are called driving flelds. Setting their F-terms to zero leads to the F-term conditions 
which control the alignment of the flavon flelds. 

The relevant effective superpotential terms of the Yukawa sector of the AF model are 

Here {■ ■ ■)r denotes the contraction of the A4 indices to the representation r. When the 
three flavon flelds ips,^,fT acquire their VEVs [T2] . 



i^s) =vs\l\ , {0 = u, (<^t) =vt\0\ , (2.2) 
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Table 1: The particle content and symmetries of the ^4 model with extra ^', ^" fiavons. 

at some high scale, the flavour structure of the Yukawa couplings is generated which, 
after electroweak symmetry breaking, gives rise to a trivial Dirac neutrino mass matrix 
mo-, a TB heavy right-handed neutrino mass matrix Mr as well as a diagonal charged 
lepton mass matrix mi. Applying the seesaw formula yields a TB light neutrino mass 
matrix which is diagonalised by f/xe = ^^tmI'*? = 0). With the charged leptons being 
already diagonal, the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix becomes 
f^PMNS = f^TB- The reason why the AF ^4 model preserves TB mixing is that the fiavon 
ips preserves the 5* generator of A^ and the absence of ^' and ^" fiavons leads to an 
accidental U symmetry, where 5, U symmetry in the neutrino sector and T symmetry in 
the charged lepton sector implies TB mixing (see e.g. [TU]). 

This result is correct at LO. It is argued in [12] that higher order operators which are 
allowed by the imposed symmetries should be considered as well, both in W^^ as well 
as in the superpotential that generates the fiavon alignments in Eq. (12. 2p . In general 
such additional terms would lead to deviations from TB mixing which are important in 
the light of the latest experimental hints for a non-zero reactor angle [HIS]. In order to 
account for a (relatively) large value of ^13 the higher order corrections need to be large 
enough which can be arranged by assuming different suppression scales for different non- 
renormalisable operators. This assumption, however, raises the question about the origin 
of the suppression scales. 

In this context, it was pointed out recently that particular ultraviolet (UV) comple- 
tions of flavour models do not necessarily give rise to all effectively allowed terms in the 
superpotential (TB]. In fact, it was found that the minimal UV completion of the AF 
model does not generate any deviations from TB mixing. The simplest way to obtain the 
desired deviations is to add new messengers to the minimal model such that certain higher 
order Yukawa operators are switched on. In the neutrino sector, the effective non-trivial 
operators allowed by the symmetries reaqj 
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{N'N'')v{^s^T)v' + {N''N%..{^s^T)v + {N''N%{^s^TK,MN'N%{i^Th , (2-3) 



all of which break TB neutrino mixing. Inspection of the flavour structure of these 
terms shows that only the first two terms lead to mass matrices that have an eigenvector 
-^(1, 1, 1)^, thus retaining TM neutrino mixing. However, requiring the messengers that 



^Note that the third term, corresponding to two independent contractions, was not listed in Eq. (32) 
of [12j despite giving non-trivial contributions to the mass matrix, i.e. contributions that differ from the 
tri-bimaximal structure. 



mediate these effective operators to be matter-like, i.e. they should have a f/(l)_R charge 
of 1, we find that we will always get a contribution from the third term of Eq. fl2.3p as well. 
Therefore, without introducing new flavon fields, we cannot find a simple UV completion 
of the AF model where the neutrino mixing features a trimaximal pattern while breaking 
the tri-bimaximal one. 

2.2 Trimaximal mixing from A4 with extra ^' and ^" flavons 

The discussion in the previous subsection leads us to consider the case where extra flavons 
^' in the 1' representation and ^" in the 1" representations of A4 are added to the model as 
already shown in Table [U This has previously been suggested (without the see-saw mech- 
anism) in \T5\ where the phenomeno logical consequences of the LO terms were studied 
numerically. However the fiavon alignment was not derived in [T5] but simply postulated. 
Remarkably, although the difference between the ^' and 1^" fiavon VEVs breaks the acci- 
dental U symmetry and thereby violates TB mixing, the presence of these flavons respects 
the S symmetry and leads to TM mixing. 

In this subsection we consider the effect on the neutrino mass matrices of adding 
flavons ^' in the 1' representation and/or ^" in the 1" representations. In the subsequent 
subsection we consider the vacuum alignment problem including these flavons. Assum- 
ing the fiavon alignments in Eq. (12. 2p . it is straightforward to find the structure of the 
charged lepton and the light neutrino mass matrices. As the charged lepton Yukawa cou- 
plings are non-renormalisable, a particular set of messengers is necessary to generate the 
required couplings. Following [TB], one can show that a minimal messenger completion 
does not generate any off-diagonal entries in mi. Therefore, the leptonic mixing matrix 
t^PMNS is solely determined by the neutrino sector. Given the symmetries of Table [H the 
corresponding renormalisable neutrino part of the superpotential is extended to, 

W^^ = yLH^N^ + {y,^s + 2/2^ + v'z^' + ylON^'N^ • (2.4) 

Inserting the fiavon vacuum alignments in Eq. (12. 2p . and assuming both ^' and ^" as 
well as SM Higgs VEVs, we obtain the Dirac neutrino mass matrix m£, as well as the 
right-handed neutrino mass matrix M/j, 



(2.5) 
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(2.6) 



where a = y,vs, /3 = y^iO, i = v'si^), l" = VsiC)- 

The complex symmetric matrix Mr is diagonalised by a unitary transformation Ur, 

UIMj^Ur = Mj^^ . (2.7) 



One can easily verify that 4^(1, 1, 1)^ is an eigenvector of MJ^Mj^, the eigenvalue being 
I (/3+7'+7") p. Therefore the matrix Ur will have a trimaximally mixed column which, for 
phenomenological reasons, will be identified as the second column. The other two columns 
are more difficult to determine and we will discuss the analytic approximate calculation 
in Section HI However, they can be parameterised as a linear combination of the first and 
third columns of the TB mixing matrix Utb- Thus Ur is of the form given in Eq. (11 ■2p . 
with the parameters ■& and p being functions of the Yukawa couplings and VEVs. Note 
that a Majorana phase matrix would still have to be multiplied on the right in order to 
render the right-handed neutrino masses real. Disregarding this and other issues about 
phase conventions, we can now rewrite the effective light neutrino mass matrix, 

mf = -m^M-^V^ = -m^U^iM'-T'Ulml = -U^ ^ f/J , (2.8) 

where the matrix Ui, is defined as 

f/. = — f/« . (2.9) 

yvu 

Due to the trivial structure of m^, see Eq. (12. 5p . it is unitary and again of trimaximal 
form. In fact, Uiy differs from Upi only in the interchange of the second and the third 
row. Assuming diagonal charged leptons, the matrix Uiy is just the lepton mixing matrix 
Uu = f/pMNS- From Eq. (12. 8p it is then apparent that the Dirac neutrino mass matrix in the 
basis of diagonal right-handed neutrinos and charged leptons, is given by m'^ = yf^f/pMNS- 
In this basis, the columns of m'^ are proportional to the columns of f/pMNS? and thus we 
see that the ^4 model with additional ^' and ^" flavons satisfies form dominance [T7] 
and hence leptogenesis vanishes [18] IQ However, form dominance does not imply that the 
effective light neutrino mass matrix is form diagonalisable, and indeed this is not the case 
since the mixing matrix Ui, depends on the Yukawa couplings and VEVs (which in turn 
determine the neutrino masses). Only in the limit where 7' = 7" do we recover a form 
diagonalisable mass matrix corresponding to the TB case. 

The size of the deviations from TB mixing are controlled by the parameters y'^ and 
1/3 as well as the VEVs {^') and {C,")- From the vacuum alignment discussion in the 
next subsections, however, there is no reason to believe that these VEVs should be small, 
i.e. one would expect a large reactor angle as well as large deviations from maximal 
atmospheric mixing. In order to meet the experimental bounds, one would therefore have 
to assume the input parameters of the flavon and Yukawa superpotential to conspire so as 
to yield only small deviations from the TB mixing pattern. Eventually this will motivate 
us to go beyond A4 to S4 where the two extra flavons ^' and C," are unifled into a doublet 
representation. 



^This conclusion can be avoided by taking into account renormalisation group corrections which effect 
the entries of the Dirac neutrino Yukawa matrix differently so that the orthogonality of the columns is 
broken 1201. 



2.3 Vacuum alignment in A4 with extra ^' and ^'' flavons 

In this subsection we scrutinise the modifications to the ahgnment mechanism of AhareUi 
and Ferugho due to the addition of fiavons ^' in the 1' and/or C," in the 1" representation 
of A4. Our main result is that contrary to [12] one need not introduce the auxihary fiavon 
singlet C, which does not acquire a VEV and so plays no direct role in the fiavour structure 
of the neutrinos. In that sense the introduction of a ^' and/or C," fiavon does not only 
serve the purpose of breaking the TB structure to a TM one, but also does not complicate 
the model at all. 

The renormalisable terms of the driving superpotential, replacing the ^ fiavon of the 
AF model by the ^' and C," fiavons, are, 

^ir" = ^T (M^T + gv^TVr) + <^5 (givs^s + 92Vsi + gWs^' + g'WsO 

H'{9i^sVs + gdi + 9^^'n- (2.10) 

As the F-term equations of y?^ remain unchanged from the AF model, we obtain - up to 
some Ai transformed solutions - the well known alignment 

(¥PT)=t^T , ^r = -^- (2-11) 

On the other hand, the F-term conditions of yj^ and ^° are slightly modified in this new 
setup. Writing {f^Si) = si and (^) = u, (^') = u' and (^") = u",we get 

fsl-S2sA /sA fsA /S2\ /0\ 

2g, si - S1S2 \+g2u[s2\ +g'su' Si + (?>" S3 = , (2.12) 

\sl-sssj \ssj VW VJ \0/ 

g4{sl + 2s2S3) + g^u^ + geu'u" = 0. (2.13) 

For simplicity, let us first consider the case where only one of the two non-trivial 
one-dimensional fiavons is present. For definiteness we assume this to be the C,'- It is 
then straightforward to work out the most general solutions to Eqs. (12. 12112. 131) . Again, 
disregarding the ambiguity caused by A4 symmetry transformations we find two possible 
non-trivial solutions, 

{^s)=vs(l], vl = -^u\ «' = -^«, (2.14) 



as well as 



3^4 g's 



Choosing the soft mass parameter mi < 0, the second solution is eliminated and the VEV 
u slides to a large scale [12]. In this way we are able to get the original alignment of ips 
and non- vanishing VEVs for ^ and ^', see Eq. (12.141) . 



We now consider the effect of having both ffavons ^' and ^" in the 1' and 1" representa- 
tions of A4. Then the terms proportional to g'^ and Qq would be switched on in the fiavon 
superpotential of Eq. fl2.1Up The corresponding extra terms in the F-term conditions of 
Eqs. f l2.12f2.13p would thus modify the physical solution of Eq. fl2.14p to 

{'^s) =vs\l\ , Vs = , u = . (2.16) 

\ij 3^4 g2 

This solution has the unpleasant feature of leading to arbitrary physics. For instance, if 
y'^u' = y'^u", where y'^ denotes the Yukawa coupling of ^" to N^N^, see Eq. (12. 4p . then this 
implies that the mass matrix Mr in Eq. (12. 6p has a tri-bimaximal structure, since 7' = 7", 
and thus the reactor angle vanishes identically. It is for this reason that models with either 
^' or C," essentially yield the same physics. However, adding both types of fiavons in A4 
generates a bothersome ambiguity in physical predictions. In the next section the above 
ambiguity is removed by unifying the fiavons ^' and C," into an 5*4 doublet r]^, whose VEV 
components are aligned along a U preserving direction, thereby restoring TB mixing, at 
least approximately. 

3 S4 model of trimaximal mixing 

3.1 The effective 5*4 model of leptons 

As pointed out in the previous section, the A4 model with ^' and/or C," fiavons cannot 
explain the smallness of the deviations from TB mixing. Furthermore, adding both non- 
trivial one-dimensional fiavons leads to an ambiguity in physical predictions. In order to 
cure these shortcomings we consider an S4 model in which the 1' and 1" representations 
of ^4 are unified into the rj^, doublet of S4 while the triplet representations remain. The 
complete list of lepton, Higgs and fiavon fields is given in Table [2]o Similar to the ^4 
model we have a t/(l)/j symmetry as well as a Z^ symmetry which separates the neutrino 
and the charged lepton sector. 

In the neutrino sector of the 5*4 model there are three fiavon fields: ipu and ^^ (analogous 
to (fs and ^ of the AF model) and rj^, (which unifies the two ^4 fiavon fields ^' and ^"). 
The neutrino part of the effective superpotential is then. 



iVfiS 



Cu 



Wsf ~ LH^N' + {i^, + i, + r^,)N'N' + ^r^,N'N\ (3.1) 



M. 



X 



analogous to Eq. (12. 4p . where, as in the A^ model, the Dirac neutrino mass matrix tak;es 
the trivial form m£, of Eq. (12.50 . However an additional fiavon Cu in the 1' representation 



^ In principle the triplet fields can either be identified with the 3 or the 3' of S'4. They differ from 
each other only in the sign of the U generator (see Appendix |^ such that all representation matrices of 
the 3 have determinant +1, while this is not the case for the 3'. In the case of the right-handed neutrinos 
A'''^ and the lepton doublet L, we are free to choose the type of S'4 triplet as long as it is the same for 
both fields, and we choose the 3. Note that the triplet fiavon ip^, must furnish a 3' of S'4 because it is 
coupled to the symmetric product N'^N'^. 
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Table 2: Lepton, Higgs and flavon fields of the Si model. 

of 5*4 serves the purpose of breaking 5*4 to A^. The U violating flavon C^y is forbidden by 5*4 
to couple to N^N'^ at renormalisable level but appears at higher order, thereby breaking 
the TB structure by a small amount, leading to a non-zero but suppressed reactor angle, 
while preserving the TM structure due to the unbroken S symmetry. Inserting the flavon 
VEVs, whose alignment is discussed later, 

(V^i.) = t;,. 1 , {riu) = ^^^ ( X ) ' ^^''^ = '^^ ' ^'^''^ = ^^ ' '^^■^) 

Wg'^ leads to TB neutrino mixing at LO, broken to TM mixing by the no n- vanishing 
VEV for the C,i, flavon at next-to- leading order (NLO). Note that the neutrino sector above 
may readily be incorporated without change into an SU{b) GUT as was done for A^ and 
5*4 models of TB mixing along the lines of [2Tti2^ . 

In order to account for the charged lepton mass hierarchy we identify the right-handed 
charged leptons with one-dimensional representations of 6*4, and distinguish them by an 
extra Z3 family symmetry, broken by a triplet flavon Lpi, as well as the doublet flavons 
rj^jfje, resulting in the effective charged lepton superpotential. 









(3.3) 



where M^. are a priori independent mass scales, although two of the mass scales will turn 
out to be equal, M^^ = M^g. Inserting the flavon VEVs, 

(^m) = ^M f 1) ' (^e) = We L j , (3.4) 

whose alignment is discussed later, and using the Clebsch-Gordan coefficients listed in 
Appendix El leads to 




TT.£.cff VnHd 

W a ~ 



^* M, 



ni 






(3.5) 



resulting in a diagonal charged lepton mass matrix m^. The larger r lepton mass is 
explained by the fact that the e and /i masses are only provided at higher order by two 
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Table 3: Driving fields and messengers of the 5*4 model. 

additional doublet flavons rj^^rje- The electron to muon mass hierarchy is accounted for by 
a hierarchy of flavon VEVs Wf. <^ w^, which can be obtained consistently assuming certain 
hierarchies in the mass parameters of the flavon superpotential. Such large hierarchies of 
VEVs are of course stable in SUSY models, and are already familiar due to the hierarchy 
between the weak scale and the flavour or GUT scales. Finally note that the charged 
lepton sector above will eventually require modification as was done for A4 and 5*4 models 
of TB mixing in order to be incorporated into an SU{5) GUT since quark mixing requires 
some small off-diagonal mixing also in the charged lepton sector [2T]-|21] . 

3.2 The renormalisable »S'4 model of leptons 

As emphasised earlier, any non-renormalisable term of an effective superpotential should 
be understood in terms of a more fundamental underlying renormalisable theory. Without 
such a UV completion of a model, higher order terms which are allowed by the symmetries 
may or may not be present. Thus a purely effective formulation would leave room for 
different physical predictions. In order to remove any such ambiguity within our S4 model 
we have constructed a fully renormalisable theory of the lepton sector. The required mes- 
sengers are listed in Table |3l together with the driving fields which control the alignment 
of the fiavons. 

With the particle content and the symmetries specified in Tables [2] and |3l we can 
replace the effective neutrino superpotential in Eq. (13.1 p by the sum of two renormalisable 
pieces, a LO piece and a messenger piece, 

Wl = IV^f ° + W^sf- , (3.6) 

where 



-"^ 



ru,LO 



W's: = yLH^N' + {y^^, + y^i, + y^'n,)N'N\ (3.7) 

^.,mess ^ xiNX,x + x\x2Vu + x,^,)N^ + M^xx' . (3.8) 

The leading order contribution to the right-handed neutrino mass matrix Mr is given 
by Wg'^ and leads to a TB structure. Given a non-vanishing VEV for the (,1, flavon, 
Wg'^'^^^ is responsible for breaking the TB structure to a TM one at higher order. The 
corresponding LO and NLO diagrams are depicted in Figure [Ho Integrating out the 



*We remark that we have suppressed the additional messenger terms ((p^ + ??j/ + ^v)x'^x'^ which are 
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Figure 1: Leading and next-to- leading order right-handed neutrino mass contributions. 



messenger pair x, x'^, the NLO diagram reproduces uniquely the effective operator in 
Eq. f l3.ip . -^rj^N'^N^. Note that a similar term with rj^, replaced by (p^, is forbidden 
by 5*4 as the symmetric product N'^N^ does not include the triplet 3. Plugging in the 
flavon VEVs and using the 5*4 Clebsch-Gordan coefficients of Appendix |A] we obtain a 
contribution to M/j of the form 
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(3.9) 



The relative minus sign arises as the C,y VEV breaks the \J symmetry of 6*4. Its presence 
leads to a deviation from the TB structure which would exist if the two matrices were 
added instead. 

The 5*4 model thus gives rise to TB neutrino mixing at LO which is broken to a TM 
mixing pattern by NLO corrections induced by the VEV of the Sa, breaking fiavon field C^v 
This naturally explains why the reactor angle as well as the deviations from maximal 
atmospheric mixing are relatively small. Provided the charged leptons are diagonal, the 
lepton mixing matrix is purely determined by the structure of the right-handed neutrino 
mass matrix M^, given in Eq. f l2.6p . with 



OL 



V\Vv 



where we have defined 



/3 = y^u^. 



7 = V?.Wv , 



7 



A 



7-A, 



7 



7 + A 



XiX2- 



WyZy 

'mZ 



(3.10) 



(3.11) 



Notice that 7' and 7" are equal at LO. The deviations from a TB mass matrix arise only 
at NLO which is parameterised by — ~ -g^. 

The charged lepton sector is formulated at the renormalisable level, using two new 



pairs of messengers, Vt^ and Vt\ 



1,2). With the particles and symmetries listed in 



allowed in W^^^^'^^^ since they become relevant only at next-to-next-to-leading order with the C^y flavon 
entering quadratically; as the Klein symmetry Zg x Zjj of the neutrino sector 19 is restored in this 
diagram by the quadratic appearance of Ci,, such a higher order term yields a TB contribution to M^. 
Therefore, the only significant term contributing to the breaking of TB to TM is the one shown in 
Figure [1] 
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H, 



^i 



H, 



^e 



Vii.Ve 



L 



—^ 



L 



i il i il i in i In 



/i' , e" 



Figure 2: S4 diagrams that lead to a diagonal charged lepton sector. 



Tables [2] and [3l we get the renormalisable superpotential for the charged leptons, 



W, 



Si 






(3.12) 



where we have suppressed all order one coupling constants. The appearance of (^y in the 
second line leads to an irrelevant correction of the Q2 messenger mass which we ignore in 
the following. Integrating out the messengers, we obtain the diagrams of Figure [2] leading 



;,eff 



uniquely to W^^ in Eq. (13. 3p . with Mq^ 



M< 



^3- 



3.3 Vacuum alignment in the S4 model 

So far we have only postulated the particular alignments of the neutrino-type flavons, given 
in Eq. (13. 2p . and the flavons of the charged lepton sector, Eq. (13.41) . In this subsection 
we first explore the driving potential relevant for the neutrino sector and prove that the 
assumed alignment can in fact be obtained in an elegant way which is similar to the ^4 
case. Then we perform a similar study of the flavon alignments for the charged leptons. 
The renormalisable superpotential involving the driving fields necessary for aligning 
the neutrino-type flavons is given as 



W. 



flavon,!^ 



96 



54 - fligifufu + g2<fuC.u + gsfuVu) + ^^{g^fufu + g^iuiu + y^i.^^ 



(3.13) 



Identifying the two components of the doublet fiavon rji, with ^" and C,', the first line 
of Eq. (I3.13P is identical to the corresponding terms in the A4 driving superpotential 
of Eq. (I2.10p . Due to the 5*4 structure, we additionally get g'^ = g'l = g^. All other 
common coupling constants gi {i = 1, ...,6) appear in exactly the same way once the S4 
and A4 indices are expanded out. As identical equations yield identical results, the first 
line of Eq. (I2.iup gives rise to an alignment which suffers from the same ambiguity as the 
alignment of Eq. (I2.16P : the ratio of the VEVs of the two doublet fiavon components, 
, remains a free parameter. This ambiguity can be removed by introducing the 



I.e. 






driving field ip^ in the 3 representation of S4. As ipi, furnishes a 3' representation, we only 
get one coupling, proportional to gj, at the renormalisable level. With ip^ already being 
aligned in the (1, 1, 1)^ direction, the resulting F-term condition enforces equal VEVs for 
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rj^^i and r]^^2- In this miniinal and elegant way we have obtained the structure of the 5*4 
ahgnments given in Eq. fl3.2p . with 






92 



ul , Wy = -^r—Uu ■ (3.14) 

25-3 



As in the AF model, the VEV u,, is driven to a non-zero value due to m? < 0. 

The remaining terms of the flavon superpotential of Eq. flS.lSp are responsible for 
driving the VEV of the (,y flavon in the 1' representation of S^. As the driving field 
D'^ is completely neutral under the imposed discrete symmetries, compare Table [21 we 
additionally encounter Higgs fields as well as fiavons of the charged lepton sector. The 
former do not play any role in aligning the fiavons as the 5*4 breaking is typically assumed 
to occur around the GUT scale. Therefore we can safely set /s = 0. This leaves us with an 
F-term condition that intertwines the d, flavon with the fiavons rj^ and rje of the charged 
lepton sector. A simple way to disentangle the fiavons of the two sectors was put forward 
in [21]. The idea is to introduce a second driving field with identical quantum numbers 
as D^. Then one obtains two F-term conditions which are identical in their structure but 
involve independent coupling constants, fi for one driving field and /,' for the other. As 
both conditions must be satisfied, one can linearly combine the equations and thus find a 
unique solution for the VEV of the d, flavon, 

4 = - Yf'Jff' ^"^ ■ ^3.15) 

/1/2 /2/1 

Turning to the flavon alignment of the charged lepton sector, the renormalisable driv- 
ing superpotential takes the form 

p^flavon,^ = iflQiiipe^ifi + h2ipeVfM) + ^eih'fe^i + ^4'7m'7m) 

+rfl{M2rie + KCuTrje + hQirj^^rj^ + hjip^ip^) . (3.16) 

The alignments of the flavon fields (pi^ and r/^ are constrained by the first line of Eq. (I3.16p . 
Writing {ip^^) = ti and {1]^.) = rj, the F-term conditions read 



+ h2n Its +h2r2 ti = , (3.17) 





h-sitl + 2t2h) + 2/i4rir2 = . (3.18) 

The most general solution to these equations can be obtained by a straightforward calcu- 
lation. As a result one finds two solutions, 

((/?^)=t;^[lj , (r/^)=w^M , ^^ = __J_y;^^ (3.19) 
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as well as 



i^e) =Vi\]] , (V,) = w, [ _M , ^' = ir ^J ■ (3.20) 



Clearly, alignments which are derived from these by application of an arbitrary 6*4 trans- 
formation also satisfy the F-term conditions of Eqs. fl3.17|3.18p ]^l In order to get rid of 
the second solution we can use the idea of adding a second ^^ driving field [21], thus 
effectively splitting the F-term condition of Eq. fl3.18p into two equations, 

{tl + 2t2h) = , nr^ = . (3.21) 

Then the remaining solution, i.e. Eq. f l3.19p . is just what we have used in the discussion 
of the charged lepton flavour structure, see Eq. fl3.4p . 

Having fixed the alignment of (fe and 77^ we can now proceed with the alignment of 
the fiavon doublet rje- This is determined from the F-term condition of the driving field 
77° in Eq. fl3.16p . Inserting the VEVs of ipi, rj^ and (^, and writing (r/g.i) = Qi we obtain 

M. (*) + hr.. (.';) + A«..J (J) + krvf (;) ^ (°) . (3.22) 

As Zi, is already fixed, 52 = and we find the rj^ alignment of Eq. ( 13. 4p with 

We = rT^ ; . (3.23) 

Finally, we have to consider the terms of W^^'^^'^ proportional to the driving field D°, 
see Eq. (I3.13p . Repeating the discussion above Eq. (I3.15p . the 77^ and ?7e VEVs are 
simultaneously driven to non-vanishing values with 



/4/1 /1/4 s/r2 
72 /l ~ /1/2 



^M^e = - fl y', K ■ (3.24) 



We conclude this section by sketching how the hierarchy of charged lepton masses 
might be obtained in our setup. We will adopt a very rough approximation, dropping 
all order one coefficients and keeping only orders of magnitude. Everything must be 
explained in terms of the input mass parameters Mq and M^ of Eqs. (I3.13f3.16p as well 
as the involved messenger masses. Assuming 

M2 ~ 1000 Mo , (3.25) 

we get 

t;, ~ w^ ~ lOMo , We ~ — Mo . (3.26) 

This yields a hierarchy of about 100 between the muon and the electron mass. Setting the 
messenger masses Mq. at a scale of around lOOMo fixes the ratio between the tau and the 
muon mass at a value of order 10. The tau mass itself is then set to be approximately j^. 



^For a discussion of why such S^ symmetry transformed solutions lead to the same physics see |23) . 
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4 Trimaximal mixing in terms of perturbative devi- 
ations to tri-bimaximal mixing 

In this section we make a perturbative and analytic study of the deviations to TB mixing 
which are predicted by the A4 and S4 models of TM mixing, subject to the phenomeno- 
logical constraint that the reactor angle is perturbatively small, which enables TM mixing 
to be viewed as a perturbative expansion around the TB limit. This is natural from the 
point of view of 5*4 models where TB mixing arises at LO, broken by higher order correc- 
tions which preserve TM mixing. It also enables alternative phenomenological proposals 
to be viewed and compared on the same footing. 

Our starting point is the right-handed neutrino mass matrix of the A4 model in 
Eq. fl2.6p . It can be written as the sum of a matrix that preserves TB mixing and a 
matrix that violates it, 

M^ = M™ + AM^, (4.1) 

where 

M™ = a\-l 2 -1|+/3|0 011+7(1 101, (4.2) 



and 




(4.3) 



with 



A = 1(7" -7'), 7 = i(7' + 7"). (4.4) 

In the 5*4 model the explicit form of AMpi is given in Eq. (13. 9p . In both the A^ and S4 
models the TB violating matrix is required to be small, 

|A| < |a|,|/3| . (4.5) 

This assumption is necessary in order to meet the experimental constraints that any 
deviations from TB mixing should be small. The parameter 7 on the other hand need 
not be small as it does not break the TB pattern. 

Since M^^ is diagonalised by the TB mixing matrix U^b, this enables Mji to be 
diagonalised perturbatively by expanding about the TB mixing case. Writing the matrix 
Uji which diagonalises the right-handed neutrino mass matrix in terms of its column 
vectors $4, 

Ur = ($l,$2,$3) , (4.6) 

we can expand Ur, for small deviations from TB mixing, in linear approximation around 
its TB form using 

$, = <I>™ + A<I>, , A$, = 5^a,,<l>™, (4.7) 
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where $™ are just the columns of the TB mixing matrix. As shown in Appendix |Bl due 
to the unitarity of Ufi and the special form of the mass matrix M/j in Eq. (14. ip . the only 
non-zero parameter is 013 = —ali whose dependence on the input parameters a,/3,7, A 
is given in Eqs. flB.15|IB.16|) . The fact that only ais = —a^i is non-zero implies that Ur 
is of TM form as expected. Furthermore, since, 



UlM^U^ 



Mf"^ 



(4. 



it is then straightforward to derive the lepton mixing matrix f/pMNS, as in Eq. (12. 9 p 

mo 



U, 



PMNS 



■u^ 



yvu 



R 



(4.9) 



Due to the trivial structure of mo as well as a diagonal charged lepton sector, the PMNS 
mixing matrix can thus be directly obtained from Upt by permuting the second and the 
third row as well as multiplying the Majorana phase matrix P on the right and another 
phase matrix P' on the left, leading to f/pMNS = f^TM where. 



t/T 



P' 



( ^ 



J_ 2_ * \ 



■775 + :75«13 775 77^ + 77S" 



V- 



V6 ' v^"^'^ yf3 V2 



1 



V2 



"13 



1 



v^"13 



P . 



(4.10) 



-71 + TS'^is/ 



The matrix P' has to be chosen such that the PMNS matrix without Majorana phases 
is brought to the standard PDG form where the 2-3 and 3-3 elements are real and the 
mixing angles are all between 0° and 90°. In linear approximation, the required form of 
P' becomes 

diag(l,a+,-a_) , a± = l±i- ^^^^' . (4.II) 



P' 



V3 



Multiplying this explicit form of the phase matrix P' we obtain a mixing matrix that is 
consistent with the standard PDG phase conventions. 

It is useful to compare the TM mixing matrix in Eq. (I4.10p to a general parametrisation 
of the PMNS mixing matrix in terms of deviations from TB mixing [25] , 



/ 



t^PMNS 



— 75(1 + s — a + re 



^^ 



\ V6 



'1 + s + a 



re 



73(1 + ^) 
''' ^(l-|3-a-|re^^) 
1 ^' ^s + a+\re'^) 



\ 



-73(1 



75(1 +«: 



P, (4.12) 



V2^ 



where the deviation parameters s,a,r are defined as [25] . 



sin 6^12 = ^(1 + s 
v3 



This comparison yields 

Re (013) 







v^ 



sin ^23 = ^(1 + « 
V2 



r cos6 



V3 



Re (a 



13 j 



sin6'i3 = —^ 



V2 



(4.13) 



5 ^ arg(ai3) + 7r , (4.14) 
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where S is the CP violating oscillation phase of the lepton sector, and 013 is proportional 
to A as shown in Eqs. (]B.15|B.16p . Independently of the value of 013 we confirm the TM 
sum rules given in j26j . 

s f» , a^ — rcos6 . (4-15) 

We emphasise that the sum rules in Eq. fl4.15p hold for any model of trimaximal mix- 
ing, since it is always the case that phenomenology requires that it has approximate 
tri-bimaximal form. The above perturbative form of TM mixing in terms of TB deviation 
parameters is useful when comparing TM mixing to other proposed forms of the mixing 
matrix. For example, if the reactor angle is measured to be sizeable, but the solar and 
atmospheric angles remain close to their tri-bimaximal values, i.e. the deviation param- 
eters in Eq. fl4.13p take the form s ~ a ~ but r 7^ 0, then the mixing matrix takes the 
"tri-bimaximal-reactor" (TBR) form [26]. Such a mixing has recently been obtained in 
an S'4 setup [2^. Alternative proposals [281 - HO] that have been put forward to accommo- 
date the T2K result could similarly be compared using the deviation parameters s,a,r. 
With future neutrino oscillation experiments being able to not only accurately measure 
the reactor angle, parametrised here as r, but also the atmospheric and solar deviation 
parameters a, s and eventually the CP violating oscillation phase S, it is clear that re- 
lating these deviation parameters via sum rules comprise the next step in discriminating 
different models of lepton masses and mixings. 

5 Conclusions 

In the well known direct models of tri-bimaximal (TB) mixing, based on A4 and 5*4, the 
TB mixing is enforced by a Klein symmetry Z|^ x Z^ in the neutrino sector, together 
with a Zj symmetry in the charged lepton sector, where a common basis corresponds to 
a diagonal charged lepton mass matrix. It is also well known that TB mixing can emerge 
from either S'4, which contains the generators S,T, U, or A4, which contains S,T. In the 
case of A4 the U symmetry emerges accidentally as a result of the absence of flavons in 
the 1' or 1" representations of A^. Such models are called "direct models" since (some 
of) the group generators remain unbroken in different sectors of the low energy effective 
theory. Although this simple and appealing picture is apparently shattered by the T2K 
results, which indicate a sizeable reactor angle 6*13, simple alternative possibilities such as 
trimaximal (TM) mixing remain. 

We have argued that, in the framework of direct A4 models, the T2K results motivate 
adding the "missing" flavons 1' or 1", whose VEVs break the accidental U symmetry 
but preserve the S symmetry in the neutrino sector corresponding to a Zf subgroup of 
A4, leading to a TM pattern of mixing as in Eq. (II. 2p . We have studied the vacuum 
alignment of such an extended Altarelli-Feruglio A4 family symmetry model including 
additional flavons in the 1' and 1" representations and show that it leads to TM mixing 
in which the second column of the lepton mixing matrix consists of the column vector 
(1, 1, 1)"^/a/3. However there are drawbacks to this approach. To begin with, the higher 
order corrections need to be kept under control otherwise they can wash out the desirable 



16 



TM structure, however this can in principle be done by formulating the theory at the 
renormalisable level. Furthermore, this approach within ^4 provides no explanation for 
why the reactor angle should be smaller than the atmospheric and solar angles. 

In order to overcome these drawbacks, we have proposed a renormalisable S4 model 
of leptons in which the 1' and 1" flavons of A4 are unified into a doublet of S4 which is 
spontaneously broken to A4 by a 1' flavon of 6*4. We have studied the vacuum alignment in 
the 5*4 model and shown that it predicts accurate TM neutrino mixing due to a preserved 
Z2 symmetry in the neutrino sector, where Zf C Z2 x Z^ C 6*4. In the S4 model, the 
U generator corresponding to Z^ is broken by the 1' flavon of S4, which however only 
enters the neutrino sector at higher order, resulting in approximate TB mixing. 

Although the A4 and S4 models of leptons presented here involve diagonal charged 
lepton mass matrices, when the models are extended to include quarks, for example in 
the framework of SU{5) unification, we would expect the charged lepton sectors (but not 
the neutrino sectors) of these models to be modified. This could introduce additional 
contributions to lepton mixing from the charged lepton sector. Interestingly both the 
A4 and 5*4 models here preserve form dominance and hence predict zero leptogenesis, 
independently of charged lepton mixing, up to renormalisation group corrections. 

It should be clear that a sizeable reactor angle as indicated by T2K does not rule 
out the family symmetry approach, on the contrary it provides additional input which 
helps to refine the symmetry approach. In particular we have explored an S4 model in 
which the leading order Klein symmetry Zf x Z^ associated with TB mixing is broken 
at higher order down to a Zf subgroup capable of enforcing a simple TM mixing pattern. 
We emphasise that the smoking gun signature of TM mixing is accurate trimaximal solar 
mixing s ~ 0, together with the sum rule 2a + rcos6 ~ 0, which can be tested in 
forthcoming high precision neutrino oscillation experiments. 
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Appendix 

A S4 and A4 Clebsch-Gordan coefficients 

Finite groups can be defined in terms of their generating elements which satisfy certain 
product rules, the so called presentation. Starting from these rules it is possible to work 
out the matrix representations of these elements up to an ambiguity related to the choice 
of basis. In the literature there exist different 'standard' choices which all have advantages 
but also disadvantages. As the Clebsch-Gordan coefficients depend on the chosen basis, 
it is necessary to carefully define the basis in which a model is constructed. 
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Table 4: The generators S, T, U of S^ and S, T of A^ as used in this article. 



The generators and Clebsch-Gordan coefficients of 5*4 = A(24) and ^4 = A(12) in a 
basis where the triplets are exphcitly real were derived in a general way in [1T1I12]. As 
was argued in [12], there exists a more suitable triplet basis for A^ models in which the 
order-three generator T is brought to a diagonal and complex form. By now this choice 
has become the standard or physical basis for direct models [19]. The corresponding basis 
in the case of 6*4 can be found for instance in [13]. As A4 is a subgroup of 5*4 it is natural 
to express this relation also in terms of the generators where the elements S, T, U generate 
5*4, while A4 is obtained by simply dropping the U generator [H]. 

Table H] lists the generators of S4 and ^4 in the physical basis. The primed representa- 
tions of 5*4 differ only in the sign of the U generator. Dropping the U generator we obtain 
A4. It is clear from the table that the doublet of 5*4 becomes a reducible representation 
under A4, denoted by 1" and 1'. 

The 5*4 product rules in the chosen basis are listed below, where we use the number 
of primes within the expression 



a 



(') 



/3W ^ 7 



(') 



(A.l) 



to classify the results. We denote this number by n, e.g. in 3 ® 3' — ?■ 3' we get n = 2. 



l(') ^ i(') ^ i(') ) n = even 




a(3 , 



iW (g) 2 



-)■ 



n = even 
n = odd 



1 

V 



2 
2 



— > 



2 
2 



a 



(3i 



lW®3W ^ 3W 



n 



even 



1 ®3 ^ 3 

V®Z' -^ 3 

1 ®3' ^ 3' 

l'®3 ^3' 




2 ® 2 ^ iW 



n = even 
n = odd 



2(8)2 -^ 1 
2(8)2 -^ 1' 



aih + (-l)"a2/3i 
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2 (g) 2 -)■ 2 < n = even 



2 (g) 3W ^ 3W < 



2^2 



"2 A 

ai/3] 



n = even 


2®3 
2® 3' 


^ 3 ^ 
^ 3' 


ra = odd 


2®3 
2® 3' 


^ 3' 
^ 3 




3(') ® 3(') ^ iW J n = even 



3W ^ 3W ^ 2 



n = even 



n = odd 



30^30 ^ 30 <^ n = odd 



3(0 3W ^ 3(0 <^ n = even 



3 ®3 

3'® 3' 
3 ®3' 




ai/3i + aa/^s + "3/32 



a2/32 + as/?! + aiP-i 
-l)"(a3/33 + ai/32 + a2f3i, 

'2aiPi - ^2/33 - a3/32\ 
2a3/33 - ai/32 - a2/3i , 
,2a2/32 - as/?! - ai/^s/ 

^a2/33 - a3/32^ 
ai/32 - a2/3i 



The y44 Clebsch-Gordan coefficients can be obtained from the above expressions by 
simply dropping all primes and identifying the components of the 6*4 doublet 2 as the 1" 
and 1' representations of ^4, see Table HI We thus find for the non-trivial products 



1'®1" -^ 1 



l'®3 ^ 3 



l"(g)3 -> 3 



a/3 




3 
3 
3 



3^1 
3^1' 
3 ^ 1" 



3®3 ^ 3+3 



ai/3i + a2/33 + a3/32 , 
as/^s + ai/32 + 012(^1 , 
a2/32 + a3/3i + ai/33 , 

^2ai/3i - 02(33 - 03/32^ 

2a3/33 - ai(32 - a2/3i 

, 2^2/32 - a3/3i - ai/33 , 



/a2/33 - ^3/32^ 
+ ai/^2 - a2/3i 
\a3(3i - ai/33 
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B Perturbative diagonalisation of Mn 

In this Appendix we present a general method for a perturbative diagonahsation of the 
right-handed neutrino mass matrix and the determination of the PMNS matrix in the di- 
agonal charged lepton mass basis. The method is analogous to the perturbative expansion 
of the effective neutrino mass matrix developed in pS] , which is mainly applicable to the 
case of a diagonal right-handed neutrino mass matrix as encountered in indirect models. 
We show that, in the case of the direct models in this paper, it is sufficient to consider 
the matrix Un which perturbatively diagonalises the right-handed neutrino mass matrix, 
since the PMNS matrix is then simply obtained from Un due to the trivial structure of 
the Dirac mass matrix. 

Writing the matrix Uji as in Eq. (14.61) . Uji = ($i,$2,$3), we can express the mass 
matrix Mr, see Eq. (14.81) . as 

Mn = U^j^mS'^'uI = J2m,^:^U (B.l) 

i 

where Mi are the eigenvalues of Mji. For small deviations from TB mixing, we can expand 
this equation in linear approximation around its TB form using 

M^ = M™ + AM^ , Mi = M™ + AMi , $i = $™ + A$i . (B.2) 

Here M^^ is the part of Mr that does not depend on the small parameter A. Its eigen- 
values are 

M™ = 3a + /3 - 7 , M2™ = /3 + 27 , M™ = 3a - /3 + 7 , (B.3) 

with eigenvectors 






Following [IS], we parameterise A$j in the basis of these TB column vectors $™, 

A$, = J]a,,$™. (B.5) 

j 

The parameters aij, which are a measure of how much the mixing deviates from the 
TB pattern, are small. In order to determine their dependence on the input parameters 
a,/3,7. A, we first observe that Ur must be unitary, i.e. 



$™ + E"^^€') (^f + E«.'^™) 



k= l$r + 2^a.,$rj l$;" + }_^a,,<l>/^l ^5.,+a^ + a,, , (B.6) 

where we have dropped second order terms. As a consequence of unitarity we thus get 

aji ^ - a*^ . (B.7) 

20 



Expanding Eq. (IB.ip around its TB structure and keeping only terms linear in small 
deviations, we get 

^ J2 (am, <l>r<l>r^ + M™ A<l>*<l>™^ + M™ •I'PA^.t) (B.8) 

i 

^ Y: ("am, $r$r^ + M™ Yl «:.*r*r^ + m™ j] $r«:,*r 



^TB-' 



With AMij corresponding to the TB breaking contribution to Mr, as in Eq. (13.91) . 

/O 1 -1\ 
AMk = a 1 -1 , (B.9) 

V-1 ij 

it is now possible to calculate the parameters aij and AMj from Eq. fIB.Sp . Counting the 
number of unknowns we get 6 real parameters for AMj as well as 3 + 6 real parameters 
for aij where we have accounted for the unitarity constraint of Eq. (]B.7|) rl On the other 
hand Eq. (]B.8|) yields 6 complex conditions, 3 from the diagonal entries and 3 from the 
off-diagonals (AM/j is symmetric). As the number of unknowns is bigger than the number 
of conditions we do not find a unique solution. Note however that the effect of an can 
always be absorbed in AMj. Hence we can remove 3 unknowns by setting an = without 
loss of generality. 

The remaining 12 unknowns can be determined by sandwiching Eq. fIB.Sp between 
$™^ and $™^ 



iT 



$™ AM^$™ ^ (AM,4/ + M™ «;, + M™ «*,) . (B.IO) 

Explicit calculation of the left-hand side gives 

^ M™a*2 - M2™ai2 , (B.ll) 

^ M2™«;3 - M3™«23 , (B.12) 

-VSA ^ M™at3 - M3™ai3 , (B.13) 

^ AMi ^ AMa ^ AM3 . (B.14) 

Notice that the TB breaking parameter A does not give rise to corrections to the mass 
eigenvalues in linear approximation. Eqs. (IB. lip and (IB.12P tell us that ai2 = 0:23 = 
unless \M^^\ = \Mj^\ and \Mj'^\ = \Mj^\, respectively. Eq. (lRl3|) can be used to 
determine the complex valued a 13 in terms of the parameters a, /3, 7 and A of the right- 



^Note that an must be purely imaginary. 
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handed neutrino mass matrix M^. A straightforward but tedious calculation yields 



Re (aia) ~ " ^ 



A \ . / A \I^ 1^ 



Re + Im 



/3-7 



(B.15) 



v^ I- A 



Im(ai3) ^ ^:^ f-^. (B.16) 

Notice that ai3 is proportional to the TB violating parameter A, which, in the 6*4 model, 
is proportional to the C,u flavon VEV as shown in Eq. (13. lip . Once ais is determined, the 
matrix Ur, and hence the PMNS matrix, are fully determined, as discussed in Section HI 
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